30th International Cosmic Ray Conference 




ICRC*/ 



Merida, Mexico 



o : 
o . 



O 

m 

si 

6 



> 

m 
in 

in 
o 
i> 
o 



x 



The inefficiency of the first-order Fermi process in UHECR production at relativistic 
shocks 

J. Niemiec 12 , M. Ostrowski 3 , M. Pohl 1 . 

1 Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, USA 
2 Instytut Fizyki Jgdrowej PAN, ul. Radzikowskiego 152, 31-342 Krakow, Poland 
3 Obserwatorium Astronomiczne, Uniwersytet Jagiellonski, ul. Orla 171, 30-244 Krakow, Poland 
niemiec @ iastate. edu 

Abstract: The question of the origin of ultra-high-energy cosmic rays at relativistic shock waves is dis- 
cussed in the light of results of recent Monte Carlo studies of the first-order Fermi particle acceleration 
[8, 9]. The models of the turbulent magnetic field near the shock considered in these simulations in- 
clude realistic features of the perturbed magnetic field structures at the shock, which allow us to study 
all the field and particle motion characteristics that are important for cosmic-ray acceleration. Our re- 
sults show that turbulent conditions near the shock, that are consistent with the shock jump conditions, 
lead to substantial modifications of the acceleration process with respect to the simplified models, that 
produce wide-range power-law energy distributions, often with the "universal" spectral index. Relativis- 
tic shocks are essentially always superluminal, and thus they preferentially generate steep particle spectra 
with cutoffs well below the maximum scattering energy, often not exceeding the energy of the compressed 
background plasma ions. Thus, cosmic-ray acceleration to very high energies at relativistic shock waves 
is inefficient, and such shocks are not expected to be the sources of ultra-high-energy particles. 



Introduction 

Relativistic shocks are widely considered to gen- 
erate energetic particle populations (cosmic rays) 
responsible for the high-energy emission of astro- 
physical sources such as hot spots in radio galax- 
ies, quasar jets and gamma-ray burst afterglows. 
The basic acceleration mechanism discussed in this 
context is the first-order Fermi process. It is be- 
lieved that the Fermi process is intrinsically effi- 
cient and thus also capable of the production of 
ultra-high-energy particles. In this work we con- 
front this opinion and show that the generation of 
UHECRs at relativistic shocks must invoke pro- 
cesses other that the first-order Fermi mechanism. 

Numerical models and results 

Modeling of first-order Fermi acceleration at rel- 
ativistic shocks is a difficult task because cosmic- 
ray distributions are highly anisotropic at the shock 
and the resulting particle spectra depend strongly 
on the essentially unknown local conditions at the 



shocks. In the series of our recent studies of the 
Fermi process [7, 8, 9] we have to considered the 
most realistic models possible for the perturbed 
magnetic field structures at the shock, which allow 
us to study all the field characteristics important for 
particle acceleration. The upstream magnetic field 
is assumed to consist of the uniform component 
Bo,i> inclined at an angle ipi to the shock normal 1 , 
and static finite-amplitude perturbations imposed 
upon it. The irregular component has either a flat 
(F(k) ~ k- 1 ) or a Kolmogorov (F(k) ~ fc~ 5/3 ) 
wave power spectrum defined in a wide wavevec- 
tor range with k max /k m i n = 10 5 , which allows 
us to investigate the role of the long-wave turbu- 
lence in the acceleration process. The downstream 
field structure is obtained as the compressed up- 
stream field, so that the magnetic field lines are 
continuous across the shock. This allows one to 
study upstream-downstream correlations in parti- 
cle motion introduced by the field structure for dif- 
ferent levels of turbulence, and to investigate the 



1. Indices "1" and "2" refer to quantities in the up- 
stream and downstream plasma rest frame, respectively. 



Cosmic-ray acceleration at relatvistic shocks 



(5B/B 0>1 = 1.0 
F(k)~k- 5/3 



7,=30 



_Lj 



7i = 5 



7-1 = 10 



12 3 4 5 6 

log(E/E ) 

Figure 1 : Accelerated particle spectra in the shock 
rest frame at oblique (?/>i = 45°) superluminal 
shock waves for different shock Lorentz factors 71. 
The Kolmogorov wave power spectrum is assumed 
for the turbulent magnetic field, and the upstream 
perturbation amplitude SB/B01 = 1.0. Particles 
in the range (2ir/k max , 27r/fc m m) can satisfy the 
resonance condition k res ~ 2n /r g {E) for some of 
the waves in the turbulence spectrum. 



influence of this factor on particle spectra. We 
study the first-order Fermi process in test-particle 
approach with the method of Monte Carlo simula- 
tions, which calculates the particle spectra by fol- 
lowing exact particle trajectories in the perturbed 
magnetic field near the shock. A shock has a pla- 
nar geometry and propagates with Lorentz factor 
71 with respect to the upstream plasma. 

Because nearly all magnetic field configurations in 
relativistic shocks lead to a perpendicular (superlu- 
minal) shock structure, the characteristic features 
of particle acceleration processes at high-7 shocks 
are best illustrated using the oblique shock exam- 
ple of Fig. 1. All injected particles are initially ac- 
celerated in a phase of "superadiabatic" compres- 
sion at the shock [2]. Only a much smaller fraction 
of these particles is further accelerated in the first- 
order Fermi process, forming an energetic tail in 
the spectrum for highly perturbed magnetic fields. 
The shape of the spectral tail and its extension to 
high particle energies strongly depend on the mag- 
netic field turbulence spectrum. The tails for the 
Kolmogorov turbulence (Fig. 1), with most power 
in long-wave perturbations, are much flatter than 
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Figure 2: Particle spectra at parallel shock waves 
with 7i=30 for different amplitudes of the mag- 
netic field perturbations SB/Bqi. Linear fits to 
the spectra are presented and values of the (phase- 
space) spectral indices a are given in italic (the en- 
ergy spectral index a = a — 2). Some spectra are 
vertically shifted for clarity. 



for the fiat wave power spectrum. However, in 
either case, the spectra steepen and/or the energy 
cut-offs occur in the resonance energy range, and 
the cut-off energy decreases with growing shock 
Lorentz factor. These spectral features result from 
the character of particle transport in the magnetic 
field downstream of the shock, where field com- 
pression produces effectively 2D turbulence, in 
which particle diffusion along the shock normal is 
strongly suppressed. In effect, advection of parti- 
cles with the downstream flow leads to high parti- 
cle escape rates, resulting in steep particle spectra. 
The existence of the Kolmogorov turbulence at the 
shock allows for the formation of more extended 
and flatter spectral components, due to the effects 
of high-amplitude long-wave magnetic field per- 
turbations which can form locally subluminal field 
configurations at the shock, thus enabling more ef- 
ficient particle-shock interactions. However, the 
importance of these effects diminishes for larger 
shock Lorentz factors. 

The effects of the turbulent field compression may 
also occur in parallel high-7 shocks (Fig. 2) for 
large-amplitude perturbations. In this case, the 
field compression leads to an effectively perpen- 
dicular shock configuration, and features analo- 
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Figure 3: Particle spectra for superluminal shocks 
with 7i=10 formed in the presence of shock- 
generated downstream turbulence and the Kol- 
mogorov power spectrum of the background 
field component (solid line — compare Fig. 1). 
The amplitudes of the short-wave perturbations, 
8B s h/ '{B 2 ), are given near the respective spectra. 



Figure 4: Accelerated particle spectra at paral- 
lel shocks with 71= 10 in the presence of shock- 
generated turbulence and weakly perturbed back- 
ground field. The spectrum indicated as ".scoff up" 
was obtained in the unphysical model in, which the 
particle pitch-angle diffusion was assumed to exist 
both downstream and upstream of the shock. 



gous to those observed in oblique shocks are re- 
covered. Only for weakly perturbed magnetic 
fields can the wide-energy range particle spectra 
be formed. However, they are non-power-law in 
the full energy range, and their power-law parts are 
flat (a < 4) due to the effects of long-wave per- 
turbations. The convergence of the spectra to the 
"universal" spectral index {a w 4.2) claimed in 
the literature [e.g., 1, 3, 4] is clearly not observed. 

More realistic microscopic models of collisionless 
shocks show that the shocks can generate a highly 
nonlinear short-wave turbulence downstream due 
to filamentation instabilities at the shock front [5, 
e.g.]. In [9] we augmented the magnetic field 
model by this shock-generated component, which 
can provide efficient particle scattering and may 
lead to a decorrelation between particle motion and 
the compressed field downstream of the shock. For 
the case of oblique shocks (Fig. 3), increasing the 
amplitude of the shock-generated turbulence leads 
to a more efficient acceleration with particle spec- 
tral tails extending to higher energies. However, 
in all cases, in which SB s h/ (-B2) ^> 1, the ener- 
getic spectral tails are convex, and the spectra have 
cutoffs at energies for which the resonance condi- 
tion for interactions with compressed turbulence is 



fulfilled. This is because the influence on particle 
trajectories of the shock-generated turbulence de- 
creases with increasing particle energy, and even- 
tually becomes smaller than the influence of the 
large-scale background field. Similar spectral ef- 
fects are also observed for parallel shocks when 
the amplitude of the long-wave background com- 
ponent is large. Extended power-law particle dis- 
tributions can be formed in parallel shocks prop- 
agating in a medium with low-amplitude of the 
long- wave turbulence (Fig. 4). However, the spec- 
tra are steeper than the expected "universal" spec- 
trum, a > a u - The only case in which we were 
able to obtain spectra with a = a u in the energy 
range considered, involved the unphysical model 
with the short-wave component introduced both 
downstream and upstream of the shock (spectrum 
indicated as "scatt up" in Fig. 4), which removed 
the effects of upstream long-wave perturbations. 

Discussion and conclusions 

Our results require a revision of many earlier dis- 
cussions of cosmic -ray acceleration up to very high 
energies in the first-order Fermi process at rela- 
tivistic shocks. The modeling shows that turbu- 
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lence consistent with the shock jump conditions 
can lead to a substantial modifications of the ac- 
celeration picture as compared to simplified mod- 
els producing wide-range power-law energy dis- 
tributions, often with the "universal" spectral in- 
dex [1, 3, 4]. The presence of highly nonlin- 
ear short-wave turbulence at the shock can lead 
to more efficient acceleration, but the amplitude 
of the shock-generated component required to pro- 
duce extended power-law spectra is unrealistically 
high, in particular for large shock Lorentz factors. 
Our simulations show that relativistic shocks, be- 
ing essentially always superluminal, possibly gen- 
erate accelerated particle distributions with cut- 
offs below either the maximum resonance energy 
enabled by the high-amplitude background turbu- 
lence (r g (E cutof f) < X(E res 

.max 

)), or approxi- 
mately at the energy of the compressed background 
plasma ions E cuto ff ~ r yim i c 2 . Thus, in conclu- 
sion, relativistic shocks are not promising sites as 
possible sources of ultra-high-energy cosmic rays. 
Should UHECR production be expected from rela- 
tivistic shocks it must invoke other processes, e.g. 
the second-order Fermi process in downstream rel- 
ativistic MHD turbulence [10]. 

Finally let us note that our models might have re- 
cently acquired observational confirmation. The 
recent Spitzer imaging of Cygnus A hotspots re- 
sulted in the detection of the high-energy tails of 
their synchrotron radiation [6]. Combined with 
data collected at other frequencies, these observa- 
tions allowed for a detailed modeling of the broad- 
band emission from the two brightest hotspots, 
which put precise constraints on the underlying 
energy spectra of ultrarelativistic electrons. The 
spectra can be approximated by a broken power- 
law with the flat low-energy spectral index ai w 
3.5 followed by a steep high-energy part with ah > 
5, with the break energy corresponding approxi- 
mately to the proton rest mass energy. Thus, the 
shape of the spectra reflects most likely two dif- 
ferent regimes of the electron acceleration process 
at mildly relativistic shocks of the hotspots: the 
preacceleration processes responsible for the spec- 
tral shape below the critical energy scale given 
by the inertia of protons, above which the first- 
order Fermi process operates. The steep slope of 
the spectra at high-energies is therefore in agree- 
ment with our modeling of the Fermi processes at 



oblique mildly relativistic shocks. In fact, the dif- 
ferences in the high-energy power-law indices and 
cut-off energies observed between the two hotspots 
may be attributed to the sensitivity of the Fermi 
process to the measured differences in the intensity 
(and possibly configuration) of the magnetic field 
at the shocks in the hotspots. 

This work was supported by MNiSW in years 2005-2008 
as a research project 1 P03D 003 29. 
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